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Abstract 

The QCD string breaking due to quark pair creation in the vac- 
uum confining field, possibly accompanied by vector, scalar or Nambu- 
Goldstone bosons, is studied nonperturbatively. The scalar light pair 
creation vertex occurs due to chiral symmetry breaking and has a 
confining form, which is computed explicitly together with subleading 
vector contributions. Dependence on light quark mass and flavor is 
specifically studied. The dominant scalar term is in good agreement 
with the 3 i"b model and experimental data. 

1 Introduction 

The topic of string breaking and effective decay Lagrangian has a long history. 
Among the nonperturbative (np) models the most famous is the so-called 
3 Po model (see pQ, [2] for details, reviews and references), which can also be 
written in the form of the interaction Hamiltonian for the Dirac quark fields 



Hi = gJ d 6 xipip (1) 

where g = 2m (? 7, m q is the constituent quark mass and 7 ~ 0.5 is the 
phenomenological parameter. In what follows we derive effective Lagrangian 
in QCD, which in the relativistic invariant form can be written similarly to 
(00) as 

C = [ ^M^d A x (2) 
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with Ai - calculated nonperturbatively through string tension, while if), ip 
are relativistic bispinors with current (polar) mass. The model ([T]) is rather 
successful (see [TJ |2j [3j H] for reviews and discussion) , however not deduced 
from the QCD Lagrangian. In the detailed analysis of heavy quarkonia in 
[5] another model was suggested, where the role of M. in fl2]) is played by 
the quark-antiquark potential Vqq(t) with confining and OGE parts of color 
octet vector type. An important step was done further in [1], where scalar and 
vector color singlet potentials were carefully studied, and similarity of results 
for 3 Po an d scalar potential model was stressed. A successful application of 
the 3 Jo model to charmonia was done in [6] with 7 = 0.4, assuming harmonic 
oscillator wave functions. This seemingly universal character of the model 
calls for the careful study of its connection with QCD. 

In the present paper we aim at the most general derivation of the string 
breaking amplitude with possible accompanying radiation of 7 or tt, 77, K, or 
vector and scalar mesons. 

Formally the problem of string breaking or of qq pair creation on the string 
connecting quark Q and antiquark Q can be reduced to the calculation of 
partition function 

Z = J DAexp(C A ) Wqq(A) det(m q + D(A)) (3) 

where La is the standard gluonic action and Wqq(A) is the (external) Wilson 
loop of (heavy) quarks QQ with fixed contour. Using world-line represen- 
tation for the det term in and retaining only one loop of quark q, one 
has 

Z lloop = J DA exp Ca (~\tr jH ^{D±z)e- K W m {A)W QQ {A)^ . (4) 

Thus the effective Lagrangian can be obtained from the averaged product 
of two Wilson loops, calculated in [7] , 

X = (W q ,W QQ ) A = ^ ex P (-aS A ). (5) 

c 

where Sa is the minimal area between contours qq and QQ. For large masses 
m q and tuq one can easily reproduce the scalar potential (sKs) model of jl] 
with £ non rci = / c"l x g — x o|<^- However in the small m q limit the transition 
from world-line representation of the light quark to the bispinor Lagrangian 
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/ ^A4ipd 4 x is not easy, since appearing of the scalar operator Ai for light 
quark implies chiral symmetry breaking. Therefore below we shall explicitly 
construct the kernel M, which exemplifies simultaneously confinement and 
spontaneous Chiral Symmetry Breaking (CSB) for quark q of vanishing mass 
m q in the field of external (possibly static for simplicity) quark Q. When 
mass m q is growing, m q 3> y/a, spontaneous CSB is replaced by the explicit 
one, and the problem is much simpler. 

To derive Effective String-Decay Lagrangian (ESDL), one can use the 
well-known Background Perturbation Theory (BPTh) [8], which was devel- 
oped further in |9j, where the confining properties of the background were 
taken into account. This theory has the prominent advantage, that it is in- 
frared safe, and the Landau ghost pole and IR renormalons are absent in the 
total perturbation theory, where the first term is purely nonperturbative and 
can be calculated via np field correlators [TUj- The latter in their turn are 
calculated selfconsistently through the only mass parameter of the theory, 
e.g. the string tension a (see [H] for a review and references). 

Therefore in computing ESDL one obtains several terms. First of all the 
well-known 3 Si term [12] due to gluon exchange, which in BPTh becomes the 
hybrid- mediated transition, discussed in [13], when confining background is 
important (for not large energy release). Below we shall be interested mostly 
in the scalar-type terms of ESDL. As was said above, those terms occur 
due to CSB in the confining background, and our derivation follows closely 
previous papers [TU dSl EH HZ] , where CSB due to confinement in QCD was 
derived in great detail. 

We shall specifically stress below the main mechanism, which creates for a 
(massless) relativistic quark simultaneously and spontaneously confinement 
and CSB, the mechanism which takes into account symmetry of quark spec- 
trum of Hamiltonian with scalar interaction. For massive quark of large mass 
the same effect occurs due to explicit CSB. Integrating out quark degrees of 
freedom in the effective Lagrangian, one obtains chiral effective Lagrangian 
for NG mesons at that point of the string, where it breaks down. In this way 
one obtains ESDL with NG meson degrees of freedom. 

The mechanism, described below and in [ISJ-P2] can be cast into explic- 
itly gauge invariant form (see appendix 1 of [16]), which is possible, since the 
systems under consideration ((Qq) or (Qq)) are white. This is in contrast to 
models, considering light quark pair (qq) alone. 

Concluding this introduction, one should stress, that our approach to 
string breaking as a CSB light-pair creation process differs in principle from 
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the Schwinger-type mechanism of pair creation, specific for vast electric fluxes 
in QED, see [IB] for discussion and references. 



2 Chiral symmetry breaking and scalar string 
breaking 

Consider now the partition function of a light quark in the field of a static 
antiquark in the presence of external currents v^, a M , s,p 

Z = J DADipDip exp [-(S + S 1 + S mt + S Q + S Q )] , 

S\ — —iJ d A xt/j^(d + m + v + 75a + s + i^pY 9 ^ 9 , 

S int = - J d*xft 'gA a t a ^ . (6) 

Here /, g are flavor indices, Sq and Sq refer to action of external quark 
currents, of (possibly high mass) quark Q and antiquark Q. 

We shall follow derivation of [HI ESI EES] , but for simplicity we shall use 
the simplest contour gauge, [IH] so-called Balitsky gauge [20J, where one can 
write 

AJx) = / aJu)FiJu)dui, a A = 1, a t = — , (7) 

JC(x) Xi 

and the contour C(x) is going from the point x = (x, x 4 ) to the point (0, x 4 ) 
on the world-line of Q and then along this world-line to X4 = —00. Note, 
that our final result (113)). ffT4]) will be cast in the gauge invariant form, which 
is the same for all contours, connecting points x, y to the world lines of Q 
(or Q). The independence of the resulting asymptotic expressions from the 
form of contours is shown in Appendix 3 of [T5] . 
Averaging over fields A^, (F^), one can write 

Z = J J D^D^exp [-{S x + S eff )} , (8) 
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where S e ff was computed in [I3]-[16j. Keeping only quadratic correlators 
and colorelectric fields for simplicitjo one obtains (for one flavor) 

S ef f = -- j d 4 xd 4 yi;(x)'j 4 [ilj(x)'ip(y)]'j i 'ilj(y)J(x,y) (9) 

where J(x, y) is expressed via vacuum correlator of colorelectric fields, [ip^] 
implies color singlet combination. Keeping only colorelectric fields, one has 

J(x, y) = ^-(A A (x)A A (y)) = f du { f d Vl D(u - v). (10) 
iv c Jo Jo 

Here D(w) is the np correlator, responsible for confinement [IT] . 

^-(F ifl {u)F kl/ {v)) = {8 ik 5^ - 5i„6, k )D(u - v) + 0{D X ) (11) 

and we have omitted the (vector) contribution of the correlator Di, contain- 
ing gluon exchange and np corrections to it. 

To start one can simplify the matter and neglect possible chiral degrees 
of freedom, contained in the 4q combination in (Q, replacing at large iV c 

Mx)${y)] -+ (^{x)^{y)) q = S q (x, y). (12) 

In this way one obtains 

S eff = J d 4 xd%ij(x)M(x,y)ij(y), (13) 

where M. is expressed via S q (x,y) 

M(x, y) = -iS q (x, y)J(x, y) (14) 

and S a is expressed via J\4, S a = - — % —r,- But one can notice, that S a contains 
the term in the denominator (m + j\4) and even for m = we look for a 



lr The fact, that quadratic (Gaussian) correlators yield dominant contribution is due to 
small vacuum correlation length A ~ 0.1 fm, which ensures small expansion parameter 
crA 2 <C 1, and is strongly supported by lattice measurements of Casimir scaling (see [TTj 
for discussion . Keeping only colorelectric correlators is justified when angular momentum 
of light quark pair qq is not large. Colormagnetic correlators produce angular momentum- 
dependent corrections \ Vcm ~ ^pr^ to the linear confinement term err, which is due 
to colorelectric correlators, see [3T] for details. 
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scalar M, which provides in S q terms with even number of 7 M , and scalar M s 
can be a self consistent solution of the nonlinear equation M s = -= — ^-^rJ ■ 

n 6 d+m+M s 

It is clear, that for m > the situation is simplified. Below we shall find 
explicitly the solution of this nonlinear equation. 

At this point one should stress, that S q (x,y) is the light quark Green's 
function in the field of antiquark. To simplify matter, we shall consider 
massless quark (m = 0) and static antiquark Q. It is clear, that without CSB 
S g contains odd number number of 7 matrices, hence trM(x, y) vanishes. To 
exemplify the appearance of CSB it is convenient to consider the limit of small 
correlation length A when M factorizes, [T4"| [15] 

M{x,y) = 6{x i -y i )M{y:,y), M(x, y) = J(x, y)/3A(x, y) (15) 

where A(x, y) is expanded in eigensolutions {ipk} of static Dirac equation 

A ( x >y) = ^2tpk(x) sign e k tp£(y), z=|x|. (16) 

J(x, y), defined in (flUj) . behaves at large and parallel |x| m |y| as 

J(x,y) w const min(x, y). 

This implies confinement, if A(x, y) ensures almost equal and parallel 
x,y. Hence for selfconsistency A(x, y) must tend to /35^{x.— y) at least for 
large |x| ~ |y|. Note, that {ipk} should be computed with the interaction 
M(x,x), which consists of the same {ipk}, as it happens in the mean field 
method. And here the spontaneous CSB reveals itself in the fact, that for 
scalar mean field M(x,x) the combination of {ijjk} entering in A(x, y) indeed 

has this property: in the 4x4 structure A = ^ ^ n ^ one obtains 

A 22 = —An (i- e - A ~ /3), if the spectrum {e n } has the symmetry property 
for e n — > —e n of the scalar potential and in this case the sum (fl6|) computed 
in the relativistic WKB method [22J yields the smeared 5^ - function at 
large x, ensuring small angle between x are y [H]. The appearance of this 
scalar structure of M ~ (3A,A ~ j3 can be called the spontaneous scalar 
generation, which gives CSB. 

In the static limit, when one can neglect the energy transfer from the 
heavy quark Q to the light quark pair qq in f JT3|) . the effective mass Ai(x, y) 
was calculated in [HI [151 [16] , using relativistic WKB approximation [22] . The 
result is 

M{x, y) = a |x| • P)(x - y)5(x 4 - y,), (17) 
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where |x|, |y| are distances to Q and the smeared 5-function 5^(z) has the 



range r = a 



x+y 

2 



while the range of 5(x^ — y±) is of the order of the 

vacuum correlation length. The explicit form of Ai(x,y) in ffT7|) is given in 
the Appendix. Hence for a long string, when r$ ^> 0.1 fm, one can neglect the 
nonlocality in ( TT3|) and take into account, that the equivalent contribution can 
be obtained, connecting the correlator J(x, y) to the worldline of antiquark 
Q, which gives the final result 



S eff = J d i xip(x)M(x)^(x) (18) 

with 

M(x) = (t(\x-xq\ + \x-xq\). (19) 
Two terms of M.{x) are shown in Figs. 1 and 2. 

There an example of trajectories (world lines) of quarks and antiquarks 
with time growing in horizontal direction. Vertical lines starting at x and 
y belong to contours C(x),C(y) used in ([7]). Cross-hatched region between 
them contributes to the kernel M(x,y) in (j!7j) . 

There are several properties of the interaction fflBl . (fT9l which should be 
mentioned. First of all, the mass term A4(x) is scalar and hence creates the 
3 Po pair qq. 

Additional terms in J\4(x), which are subleading at large distances, were 
studied in [HI [151 EE] see Appendix below. In particular, the contribution of 
the term Di in (fTTI) . which was neglected above, yields terms proportional to 
(3 in Ai(x), which are of the relative order O (jjtjk) , r — |x — xq|, and hence 
can be neglected at large r. 

Secondly, the kernel f|T9|) does not depend on flavor, which can be checked 
experimentally, as discussed in concluding section. However, the local form 
( fl9l) is an approximation of the nonlocal one, given in ( fl7|) and in Appendix, 
Eq. ( 1A.7I) derived for zero mass m q . The effect of nonlocality is numerical 
reduction of M.(x) ( Tl9l) by (10-=-15)%, as can be seen from Fig.l of [T5] . 
For larger m q the range of nonlocality decreases, as shown in Appendix, Eq. 
( lA.lOj) . and M(x, y) tends to the final form ([191 . 



The potential kernel ( !T9|) has similarity with the kernel, suggested by 
Eichten et al [5] , however, there was assumed its color octet and vector char- 
acter, i.e. A4(x) was taken to be proportional to 74^, while in our case it is 
scalar and color singlet. This assignment of [5] is not favored phenomeno log- 
ically, see [3J H] for discussion. The kernel f fl9|) coincides formally with the 
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so-called sKs kernel, which was studied nonrelativistically in [6], and shown 
to be phenomenologically successful and close to the 3 P model. 



3 Emission of mesons, accompanying the string 
breaking process 

We return here to the 4q effecting action (j9]), and will follow the procedure of 
[IT] , where chiral Lagrangian was derived in conjunction with the confining 
kernel (TTTj) . Doing bosonization in (Q with bosonic scalar variable M s and 
pseudoscalar <p a , one obtains the effective quark-meson Lagrangian 



S, 



Z = J DipD4)DM s D(j) a exp [Sqm] , (20) 
om - i d 4 xd 4 y ^ f a (x)(i(d + v + 75a + s + i^ 5 p)^8^\x - y)+ 

+ tM s (x, y)Ul${x, y))tf(y) - 2N f {J{x, y))" 1 M*(x, y)} , (21) 

U%(x, y) = exp {ij 5 t n cf> n (x, y))% . (22) 

Integrating out quark fields, one obtains effective chiral Lagrangian given 
in [TB] for the field (f) a with external currents included. 

Finally, classical equations of motion define the stationary point condi- 
tions, 

^\x,y)=0, (23) 



Mf\x,y) 



-J(x,y)Tr f4 (S q (x,y)) , 



(24) 



S q (x,y) = SV(x,y)| 0=o . (25) 

where Trj d is the trace over flavor and bispinor indices, S$(x,y) is the total 
quark propagator with chiral and extra mesons included, 

1 



S<p(x,y) = ( x 



/J 



and 



iA + iM s ^ ta ^ a 
A = d + m + v + 75S + s + ijsp 



(26) 
(27) 
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Note, that the equation (124j) coincides with f ll4p . when in S external 
currents are retained. Several properties of the new effective action §2T§ are 
to be noted. 

First of all, chiral degrees of freedom are now contained in the new string 
breaking term, which in the local approximation has the form, 

S^f = -if *xfc(x)M(x)zm%(x)tf(x). (28) 

Here /, g are flavor indices, a, (3 are 4-spinor indices, and color indices of 
ijj and ip are suppressed, while M(x)U are color-blind. The factor in 
( 128]) takes into account both np and perturbative renormalization of the NG 
fields, which are nonlocal and free (without interaction) in (|2T|) . ( 122|) . while 
in (f2"8"|) already physical local fields are contained. 



U = exp(z7 5 a t a ) = exp ( 275^^ ] , (p a X a = \Pl 



( JL + n + k+ \ 



- JL - 2EL K° 

" ' V6 v/2' n 



(29) 

In a similar way the string breaking with emission or absoption of a non 
- NG meson is described by the first term in ([2"TT) . see Fig. 2 

= -»' / d^lixYfyixWpix), (30) 

where 

5*5 (^) = (« + 7sO + s + «75P)5- ( 31 ) 

At this point one should stress that the decay of QQ meson into (Qq) (Qq) 
mesons is possible only due to string breaking, while decay into (Qq)(Qq)+ 
light meson is going via the string breaking emission mechanism of Eqs. (|28|) . 
( |30|) and in addition via two-step process with string breaking and subsequent 
emission from light quark (antiquark) in heavy-light products. When Q(Q) 
is a light quark, there is in addition a possibility of emission of NG meson 
(or light vector meson etc.) directly from Q. 
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Fig. 2 



4 Matrix elements of string decay and hadron 



emission 



In this section we shall study physical matrix elements based on effective 
actions (115) , p8]) and (j3"D]). 

The basic quantity, which can be calculated from the effective Lagrangians, 
derived above, is the two-step amplitude, which incorporates string breaking 
transition to intermediate two-body states and back to QQ state. 

In this way the transition amplitude from the state (QQ) n to (QQ) m via 
string breaking to the intermediate states (Qq)n 2 (Qq)n 3 with energy E naTl3 can 
be written as (see [231 I2H 125] for details, a similar amplitude was introduced 



■w, 



7 J+ 

° nn 2 nz° ran 2 nz 

(2n)\^ s E-E n2n3 (p) 



(32) 



Here J„ 



J, 



p) is the overlap matrix element with M(x) — > M(q) 
d 3 q 



nn 2 n 3 



J/123 7^#+(P + q)M(q)^ n2 (q)^ n3 (q) 



(33) 



and ?/i23 is the trace of normalized spin-tensors corresponding to spin-angular 
parts of meson states, while \l/ n , ip n are the radial parts, tables of I/123 are given 
in [231 125], e.g. for the 1 state n one has 7/123 = ^(q% — 2(u+n) ) wnere w > ^ 
are average kinetic energies of light and heavy quark respectively, tables of 
u, Q are given in Appendix 1 of |23j . 
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For the string breaking operator fll9p the overlap integral (1331) has the 
form 



J, 



nn 2 n 3 



JP 



O / 2/123 



"' q ^(p + q) 



(2?r 



dq 



#n 3 (q) 



Vv>(q) 



(34) 

Since heavy-light meson wave functions are well reproduced by a single 
oscillator function [231 1241125]. 



^n 2 (q) 



3/2 _ q 2 



e 2 ^ 



the resulting effective vertex M(q) can be written as 

AUq) = 2^ « £ S M, 

P2 P2 

This constant parameter M w was systematically used in 
calculate decays of bottomonium. In this case ^(B) = 0.48 GeV and m 
0.8GeV; this value was exactly used in [231 I2U 125], and one can see, that it 
agrees with the vertex mass operator M(x). 

Using w nm , one can calculate production and decay of all states involved. 



(35) 
to 



5 Comparison to 3 Po an d jKj models 

Our derivation of the effective Lagrangian f|T9|) implies, that the basic inter- 
action behind the string breaking is the scalar confining term acting between 
the light and heavy (anti) quark, (with subdominant terms of vector color 
singlet OGE type etc.) This is in exact correspondence with the so-called 
sKs model, studied in detail in [6] and compared to the 3Po model, where it 
was shown that both are successful and give similar results. In this way our 
derivation of M(|x|) in (fl~9|) gives additional theoretical status for both sKs 
and 3 Po models as candidates for the string breaking kernel. The difference 
between our approach and sKs model in [B] lies in the treatment of quark 
motion in the decay matrix elements, which was taken nonrelativistic in [6], 
with constituent quark mass m q , while in our approach we are using relativis- 
tic formalism for light quarks with current (zero) mass. In particular, our 
formalism for calculation of the factor ?/ 12 3 via vertex Zi factors in Appen- 
dices 1 and 2 of [22] and the use of relativistic string Hamiltonian for (QQ) n 
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and {Qq) n states, takes into account relativistic kinematics. With all that 
difference, results of our analysis are qualitatively and even quantitatively 
similar to those of sKs and 3 Po models. 

In particular, the width and shift of the state (QQ) n is given by 

(AE n , T n ) = (Rew nnj -Imw nn ) (36) 

and for the decay T(4S') — > BB one obtains 

T AS {BB) = 0.0033 |<WP)| 2 (GeV). (37) 

One can now estimate M w from the width of the decay T(AS) — > BB, as it 
was done in [24J. The wave functions of all participants have been found with 
good accuracy from the relativistic string Hamiltonian [26] and parametrised 
by 15 oscillator functions. Comparing experimental value r ex p(T(4S') — > 
BB) = (20.5 ±2.5) MeV [27] with flU} one obtains M w = 0.8 GeV. Another 
check of our kernel (fl9|) and its average value M w is the decay of -^(3770) 
into DD with the width T oxp = 25.4 MeV [27]. Following [25], one takes this 
state as l 3 Di and approximates it with 5 oscillator functions. Calculating 
(]33j) with M{q) M w = 0.8 GeV, one obtains Ija = 22 MeV, which agrees 
with r exp within experimental (~ 10%) accurac}o 

One can now compare this value with ffl~9]) and take into account, that 
the r.m.s. radii for T(4S) and B are 0.9 fm and 0.5 fm [26] respectively. 
Therefore, one obtains airs ~ 0.9 GeV, and art ~ 0.81 GeV, in good 
agreement with (|35|) . 

In a similar way the decay of T(5S') into six channels of BiBk, where 
Bi(Bk) = B, B* , B s , B* was considered in [28J. The total width computed in 
[2H] is r tot = 116 MeV Ml, where M u is in GeV. Comparing to r tot (exp) = 
(110±13) MeV [27], one has M u = (0.91-^1.03) GeV in reasonable agreement 
with (13 5p . Also the experimental ratio of decay into beauty-strange mesons 
is well reproduced in [28], which supports the flavor-blind kernel M in ( jT4"|) . 

We now turn to the string-breaking emission process. It is best studied in 
the transitions of the type (QQ) n — > (QQ) n , ' n '7r [231 12H [28] and in (QQ) n — > 

In [221 121] the subthreshold string breaking was considered for the decay 
of T(4S'), while in |28j the decay of T(5S") was analyzed, In all cases of this 

2 The author is grateful to V.D.Orlovsky, who provided detailed calculations of this 
decay (to be published) . 
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double-pion or single eta (n, n') transitions the effective factor appears 
to be small, Z« « « O (^). 

In conclusion, we have derived the relativistic string breaking kernel, 
which is the scalar color-singlet confining interaction, which is flavor blind 
and nonlocal for zero mass qq pair, tending to the confining or potential for 
long breaking string. This form is close to 3 P and sKs models in the nonrel- 
ativistic formalism. The author is grateful to D.V.Antonov for many useful 
remarks and suggestions. Financial support of RFBR grant no. 09-02-00 620a 
is gratefully acknowledged. 
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Appendix 

Properties of the mass kernel, Eq.(17) 

To simplify calculations one can choose the correlator D(z) in the Gaus- 
sian form, D(z) = D(0) exp (— j^j , where A is the vacuum correlation 
length, A « (0.1 0.2) fm, then a = \j D(z)d 2 z = 2ttD(0)A 2 . In the fi- 
nal expressions for the asymptotics of .M(x, y) only a enters, therefore the 
value of A and the form of D(z) influences only the small x,y region (the 
same is true for the dependence on the contours of the contour gauge, see 
Appendix 3 of PS]). Now the mass operator of Eq. ( fT5l) has the form 

M(p A = 0,x,y) = M(x,y) = a(xy)/(x, y)/9A(x, y), (A.l) 

where 

1 r 1 r 1 1 
/( x >y) = 7T7=7 / ds / dtexp[-(xs -yt) 2 },x,y = —(x,y) (A.2) 
ZWttA Jo Jo IX 
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with asymptotics /(x, x) = •i, |x| — > oo, and 

A ( x >y) = ^i(x)siffri£^+(y), (A.3) 

n 

while ?/> n (z) satisfy equations 

(t^ + /3 / +/3 1 - M ( z ' w )^n(w)rfw = e n ^n(z). (A.4) 



The mean- field-type equations flA.ll) . (IA.4[) have been solved in [Tl]~[i~6]. 



using relativistic WKB method [22], and below are listed expressions for the 
asymptotics of M(x, y), extracted from [IB], M(x, y) is a 4 x 4 matrix, which 
can be represented in the 2x2 form, with entries expressed via Pauli matrices 

« a I \ ( -Mil M-Yl \ . . 1 , / » _\ 

In the region where cos 9 between x, y is close to one, and |x| ~ |y|, one 
obtains 

a n = <t|x|Ai(z, y, 9),b u = --^-^x, y, 9) (A.6) 

47T 

2 2 

a 12 = £-Ai 2 , bia = ^-(nA' 12 + (n x L)A? 2 . 

The symmetry of M ik is : M 2 2 = M u ,M 2 i = —M\2- 
Here (|x| = x, |y| = y) 



a 2 x 2 KAax\ (x — y) 2 + 6 2 x 2 ) 

A 1 (x,y,9) = U , Vl (A.7) 

2n 2 J{x - y) 2 + 9 2 x 2 



A[(x, y, 6) = K \axy^9 2 + • (A.8) 

Note, that Ai plays the role of smeared normalized function <r 3 )(x — y) : 

J A 1 (x,y,9)d 3 y = l, (A.9) 
where one is using relation f 5(1 — cos 9) d cos 9 = ~. 
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All functions A 12 , A' 12 , A'^ are antisymmetric in x,y and therefore are 
zero in the local limit x = y, therefore they are not given here, see [To] . 
Expressions (1A.7|) . (I A. 81) are calculated for the case m q = 0, however 



the derivation in [TU [16] is valid also in the case m ^ 0, where one can 
approximately at large x replace ax in flA.7j) by ax + m q , so that Ai becomes 

(ax + m ) 2 Ki((ax + m a )\ (x — y) 2 + 8 2 x 2 ) 
Ai(x, y, 6) -> qJyK =J+ -. (A.10) 

27i 2 ^J(x- y y + e 2 x 2 



One can see in flA.lOj) . that the range of nonlocality, \x — y\ e ff, i n the 



limit of large m g tends to zero, and one can replace Ai by <5( 3 )(x — y) This 
leads to a moderate increase of effective confinement with growing m q . The 
same sort of mass dependence of confining interaction between quark and 
antiquark was observed recently on the lattice [3D] . 
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